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I/"") , Abstract. We present the results of a monitoring campaign of three eclipsing 

high-mass X-ray binaries (HMXBs: SMC X-l, LMC X-4 and Cen X-3). High- 
resolution VLT/UVES spectra are used to measure the radial velocities of these 
systems with high accuracy. We show that the subsequent mass determination 
of the neutron stars in these systems is significantly improved and discuss the 
O ' implications of this result. 

CO 

> : 

■ 1. Introduction 

>< 

A neutron star is the compact remnant of a massive star (M^ 8 M©) with a 
central density of 5 to 10 times the density of an atomic nucleus. The global 
structure of a neutron star depends on the equation of state (EOS) under these 
extreme condition s, i.e. the relation between pressure and density in the neu- 
tron star interior J Lattimer fc Prakash lEoolh Given an EOS, the equations of 



hydrostatic equilibrium can be solved resulting in a mass-radius relation for the 
neutron star and a corresponding maximum neutron-star mass. The "hardness" 
of the EOS depends e.g. on how many bosons are present in matter of such a 
high density. As they do not contribute to the fermi pressure, their presence 
(e.g. iBrown fc Bethel 119941 ) will tend to "soften" the EOS. For a soft EOS, one 
of the astrophysical implications would be that n eutron stars cannot hav e a large 
mass (e.g. < 1.55 M for the EOS applied bv IBrown fc Bethel Jl99l)); for a 
higher mass, the object would collapse into a black hole. 

Given the above, the accurate measurement of neutron-star masses is es- 
sential for our understanding of the equation of state (EOS) of matter at supra- 
nuclear densities. Just by finding one neutron star with a mass higher than the 
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maximum mass allowed by a given EOS proves the invalidity of th at EOS. Cur- 
rently the most massive neu tron star is the X-ray pulsar Vela X— 1 ( Barziv et al.l 
120011 : lOuaintrell et al.ll2003T ) with a mass of 1.86 ± 0.16 M . The milli second 
radio pulsar J0751+1807 may have an even higher ma ss: 2.2 ± 0.2 M ()Stairsl 
l2004h . Both results are in favor of a stiff EOS (see also ISrinivasa 

Another issue is the neutron-star mass distribution: the detailed supernova 
mechanism producing the neutron star is not understood, but it is likely that the 
many neutrinos that are produced during the formation of the ( proto-) neutro n 
star in the center of the collapsing star play an important role (|Burrowsl l2000). 
iTimmes et al. (1996) present model calculations from which they conclude that 
Type II supernovae (massive, single stars) give a bimodal neutron-star mass 
distribution, with peaks at 1.28 and 1.73 M Q , while Type lb supernovae (such as 
produced by stars in binaries, which are stripped of their envelopes) will produce 
neutron stars within a small range around 1.32 M . The massive neutron star 
in Vela X-l would be consistent with the second peak in this distribution. 

Neutron stars are detected either as radio pulsars, single or in a binary 
with a white dwarf or neutron star companion, or as X-ray sources in binaries 
with a (normal) low-mass companion star (LMXB) or a high-mass companion 
(HMXB). Most HMXBs are Be/X-ray binaries, but we focus on the initially 
more massive systems: high-mass X-ray binaries (HMXBs) consisting of a mas- 
sive OB supergiant and a neutron star or a black hole (Kaper & Van der Meer, 
these proceedings). In five of these systems containing an eclipsing X-ray pulsar 
the mass of the neutron star has been determined. The masses of all but one 
(Vela X— 1) are consistent with being equal to 1.4 M . However, most spectro- 
scopic observations used for these mass determinations were carried out more 
than 20 years ago, before the advent of sensitive CCD detectors and 8m-class 
telescopes; therefor e, the uncertainties are to o large to measure a significant 
mass difference ( see IVan Kerkwiik et al. lEH). 



2. Observations 

Here we present new, more accurate determinations of the mass of the neutron 
star in three of these systems, i.e. SMC X— 1, LMC X— 4 and Cen X— 3 using 
the high-resolution Ultraviolet and Visual Echelle Spectrograph UVES on the 
Very Large Telescope (VLT). These systems are in a phase of Roche- lobe over- 
flow, have well determined, circular orbits (P or b of a few days), and an optical 
counterpart of V — 14 mag, i.e. well within reach of VLT/UVES. 

Cen X— 3 is located in ou r own galaxy and was identified as the first X- 
ray pulsar IChodil et alJ (Il967f) : iGiacconi et aD ( 197lT ) . The system consists of 



an 06-7 II-IH star (|Ash et al.l Il999h with a neutron star in a 2.09 day orbit 
( Nagase et alll992ft . SMC X-l is locate d in the Small Magellan ic Cloud (SMC). 
The system hosts a B0. 5 lb supergiant (IHutchings et al.lll97"^ l with a neutron 
star in a 3.89 day orbit ((Nagase et al.lll992l ). LMC X— 4 is located in the Large 
Magellanic Cloud (LMC) and consists of a ~ 20 M^ 08 III sta r (Kaper et al., 
2005) and a neutron star in a 1.41 day orbit (|Levine et al.ll2000h . 

We have obtained a dozen spectra of all three systems with VLT/UVES in 
service mode in the period October 2001 to March 2002. The spectra have a 
resolving power of ~ 40 000, sufficient to measure line positions with an accuracy 
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Figure 1. Spectrum of SMC X— 1 in the wavelength range 3750-4250 A, 
obtained by the blue arm of the UVES spectrograph at orbital phase <f> ~ 0.51. 
Many H and He lines are detected as well as a few metal lines. The spectrum 
hasaS/N-45. 

of a few km s" 1 and cover the wavelength range 3600-6600 A. Fi gure n shows 
part of the spectrum of SMC X— 1. The spectra of LMC X— 4 and Cen X— 3 are 
of similar quality. For an extensive description of the full dataset, we refer to 
Van der Meer et al. (2005). 

3. Spectral Analysis 

To obtain a radial velocity measurement, often the complete spectra are cross 
correlated with a template spectrum. In our case the spectra are of such high 
quality that a radial velocity can be determined for each line separately. The 
advantage of such a strategy is that it is possible to map the influence of possible 
distortions due to X-ray heating, gravitational darkening and geometry of the 
line formation region. Here we select three lines, i.e. H I at 3797.90 A, He I at 
4026.21 A and He I at 4471.50 A, that are not affected by the stellar wind and 
that can be clearly separated from other lines. In Van der Meer et al. (2005) 
we present a velocity moment analysis to select the lines with objective criteria. 
To determine the line centre we fit the lines with a gaussian profile. In this way 
we obtain a radial velocity curve, which we fit with a sinusoidal profile (i.e. a 
circular orbit). 

The X-ray pulsar's orbital period is accurately known from pulse time de- 
lay measure ments. To determ i ne the orbital phase of the syste m we use the 
ephem eris of IWoidowski et alJ (jl99ST > , iLevine et alJ ((20001 ) and iNaease et al.l 
(1992) for SMC X-l, LMC X-4 and Cen X-3, respectively. The orbits are 
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circular. The fit parameters that remain are the amplitude of the radial velocity 
curve, Ko, and the 7 velocity of the system, t> 7 . A few example radial velocity 
curves are shown in Fig. [21 for the H 1 line at 3797.90 A and the He 1 line at 
4471.50 A. 



4. Neutron Star Mass 

In order to measure the mass of the neutron star and its optical companion we 
apply the mass function. For a circular orbit it can be shown that this is defined 
as: 

where Mo and Mx are the mass of the optical component and the X-ray source, 
respectively, Ko and Kx are the semi-amplitude of the radial velocity curve, 
P is the period of the orbit and i is the inclination of the orbital plane to the 
line of sight. The value for Kx and P can be obtained very accurately from 
X-ray pulse timing delay measurements. Our observations provide a value for 
K g. For the determination of the inclination of the system we use the approach 
of Rappaport k, Jos'sl (fl9^ 83'l. who showed that: 
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sini« (2) 

cos v e 

In this form the equation only holds for an optical component filling its Roche 
lobe, where is its Roche-lobe radius, a' is the separation of the centres of 
masses of the two components, and 9 e is the semi-eclipse angle of the compact 
object. Since in most of these systems soft X-rays are absorbed by the extended 
stellar w ind of the optical c ompanion, the eclipse las ts longer at low energies 
(see e.g. lHaberl et al 1 119941 : IVan der Meer et al. 1 12004ft . A more accurate value 
is obtained from high energy measurements. The ratio of the Roche-lobe radius 
and the orbital separation can be approximated by: 

— £ w 0.376 - 0.227 log - 0.028 log 2 — ^ (3) 
a' Kx Kx 

The value of the constants were determined bv lRappaport &: Josd (|l983t ). These 
depend on the ratio of the rotational period of the optical companion and the 
orbital period of the system. For a system undergoing Roche-lobe overflow, i.e. 
the system is tidally locked, one can assume that this ratio equals 1. This can be 
verified by determining the rotational period of the optical companion separately 
by measuring its (projected) rotational velocity. 

Using the formulas above we determine the mass of the neutron stars in our 
sample. The results and literature values are listed in Tabled 



5. Conclusions 



We show that with our VLT /UVES observations the mass measurements of the 
neutron stars in HMXBs are significantly improved. The masses are not all 
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Figure 2. For each selected line that is clearly separated, that has a unique 
identification and is not affected by the stellar wind of the OB-supergiant a 
radial velocity can be determined. Here the velocity amplitudes of the H I line 
at 3797.90 A and the He I line at 4471.50 A with their fitted radial velocity 
curve is shown for SMC X— 1 and LMC X— 4. The lower panels show the 
residuals of the fit to the datapoints. 



higher than 1.4 Mq, as is the case for Vela X— 1. The mass of SMC X— 1, which 
is 1.05 ± 0.09 Mq, is actually the lowest neutron star mass measured so far. 
Therefore, we conclude that a mass distribution of neutron stars is present in 
HMXBs. Possible distortions by gravitational darkening and X-ray heating, that 
can influence the value of the radial velocity amplitude of the optical companion, 
are discussed in Van der Meer et al. (2005). 

The most accurate neutron star masses have been derived for the binary 
radio pulsars. iThorsett h Chakrabartvl (^99) showed that most of these pulsars 
have a mass that is consistent with a small range near 1.35 M . Currently, 
the sample contains m ore systems a nd it is clear th a t a b roader range from 
1.25 - 1.44 M exists (|Stairsll2004h . I Van den Heuvell (|2004h suggests that the 
secondary formed neutron star in these systems tends to be the less massive of 
the two. On the other hand, our result shows that the primary formed neutron 
star can have a low mass as well. A mass distribution of neutron stars in binary 
radio pulsars as well as in HMXBs sets an important constraint on the formation 
mechanism of neutron stars. 
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Table 1. 


List of all the parameters obtained from literature. 


All errors are 1 a. 


system 


P (days) 


a x sin i (It 


s) e c 


K x (km s- 1 ) 


SMC X-1 
LMC X-4 
Cen X-3 


3.89229090(43) 
1.40839776(26) 
2.08713845(5) 


53.4876(4) 
26.343(16) 
39.56(7) 


29 ±2 
27 ±3 
34 ±3 


299.595±0.002 
407.8 ±0.3 
413.2 ±0.7 




system 


K (km s- 1 ) 




M (Mo) 


Mx (M ) 


SMC X-1 
LMC X-4 
Cen X-3 


20.3 ±0.9 
34.3 ±1.2 
26.0 ±2.5 


68 ±3 
65 ±4 
73 ±10 


15.5 ± 1.5 

15.6 ± 1.8 

19.7 ±4.3 


1.05 ±0.09 
1.31 ±0.14 
1.24 ±0.24 
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